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SUMMARY 


During  May  1987  the  author  presented  a  paper  "Range-Pair  Exceedances  in 
Stationary  Gaussian  Processes"  to  the  First  DIP  Conference  "Reliability  and 
Optimisation  of  Structural  Systems”  at  Aalborg,  Denmark. 

This  memo  describes  this  Conference,  the  better  known  ICASP5  and 
discussions  at  six  establishments  visited  during  the  same  trip. 
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1.  INTRODUCTION 

In  1986  the  author  was  invited  to  submit  a  paper  to  the  "1st  IFIP 
Working  Conference"  on  "Reliability  and  Optimisation  of  Structural  Systems" 
[11  at  the  Institute  of  Building  Technology  and  Research  at  the  University  of 
Aalborg,  Denmark.  DTP  is  the  acronym  for  the  International  Federation  of 
Information  Processing  and  the  Conference,  between  6th  and  8th  May,  1987, 
was  organised  by  their  technical  committee  TC-7  chaired  by  Professor  P. 
Thof  t-Christensen . 

Visits  were  also  arranged  to  six  fatigue  establishments  to  discuss 
work  in  fatigue,  especially  fracture  mechanics  and  crack  growth.  By  a 
fortunate  opportunity  I  was  also  able  to  attend  the  Fifth  International 
Conference  on  Applications  of  Statistics  and  Probability  to  Soil  and  Structural 
Engineering  (ICASP5).  Many  delegates  attended  both  conferences.  The 
subject  matter  (21  also  overlapped  but  visiting  both  of  them  reinforced  this 
while  presenting  a  broader  view  of  structural  reliability. 

1.1  Background 

During  1985/86  I  have  been  developing  a  program  M-FREL  for  computing 
fatigue  life  distributions  of  structures  with  several  cracks  which  interact  in 
non-trivial  cases.  This  is  based  on  the  theory  outlined  at  the  Weibull 
Conference  (3,41  but  computer  capabilities  have  allowed  significant 
extension.  The  executive,  database  and  time-stepping  parts  of  M-FREL  are 
virtually  complete  and  incorporate  standard  two-stage  fatigue,  damage 
tolerance  with  fixed  and  randomised  initial  cracks  and  intermediate  models 
allowed  by  the  presence  of  several  cracks  [51. 

The  remaining  lower  level  segments  to  be  created  include  a  menu  of 
cracking  and  damage  laws,  details  of  crack  interaction  procedures  and 
attrition  or  risk  of  overload  fracture.  The  invitation  to  the  IFIP  Conference 
was  therefore  an  excellent  opportunity  to  discuss  particular  cases  at  this  and 
at  other  establish-ments.  Because  it  was  apparently  further  developed  I 
decided  to  present  some  results  about  range-pairs  at  the  DTP  Conference  [61. 
At  ICASP5  there  was  more  interest  in  range-pairs  [71  although  I  had  some 
discussion  with  Steen  Krenk  at  Aalborg. 
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2.1  Institute  of  Aeronautics,  University  of  Pisa  (Friday,  24th  April,  1987) 

My  main  purpose  here  was  to  discuss  the  results  of  their  work  (Grant 
DA  ERO-78-G-107)  for  the  United  States  Army  and  similar  more  recent  work 
on  fatigue,  fracture  mechanics  and  crack  growth  for  built-up  stringer-sheet 
structures.  This  has  been  discontinued  but  discussions  with  Professor  Salvetti, 
Dr’s  Lazzeri  and  Frediani  greatly  clarified  their  reports  18,91.  Professor 
Salvetti  said  that  this  would  be  updated. 

These  discussions  were  arranged  during  a  preliminary  visit  which  also 
scheduled  some  time  for  the  ARL  reliability  and  range-pair  work.  On  the  day 
originally  planned  Salvetti  was  at  Aeritalia  in  Naples.  The  Institute  is  acting 
as  a  consultant  for  the  AMX  and  ATR42  projects  which  are  also  associated 
with  Aeromacchi  and  Embraer  of  Brazil. 

Their  particular  interest  is  in  the  design  of  bonded  panels.  On  the 
Wednesday  I  was  shown  some  of  the  testing  for  these,  other  industrial 
contracts  and  some  aluminium-lithium  specimens. 

The  University  work  discussed  on  Friday  24th  was  mainly  about 
cracking  of  stringer  sheet  structures.  There  was  also  a  short  description  of 
Frediani’s  [101  use  of  isoparametric  co-ordinates  for  computing  3-D 
J-integrals  by  finite  element  methods.  The  procedure  is  elegant  and  robust 
with  respect  to  element  size.  It  is  known  that  along  a  crack  front  the  J- 
integral  is  the  sum  of  a  surface  integral  and  a  local  Rice  type  of  line 
integral.  When  these  are  computed  in  terms  of  co-ordinates  for  20-node 
elements  it  is  found  that  1/4-point  or  extremely  small  elements  are  not 
needed.  Crack  fronts  are  placed  at  the  joins  of  4  adjacent  elements. 

Poe’s  fill  original  model  for  stress  intensities  in  reinforced  sheets 
was  restricted  to  concentrated  stringers  with  no  allowance  for  out-of-plane 
bending.  Work  at  Pisa  began  in  1973  with  finite  element  analyses  [81.  This 
was  followed  by  extensions  of  Poe’s  model  to  incorporate  2  bending  moments 
in  stringer  stiffnesses  for  an  effective  area  181.  The  Douglas  model  for 
stringer  flexibility  was  also  improved  and  it  was  discovered  that  this  was 
significantly  affected  by  the  elongation  of  holes  during  fatigue.  In  an 
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application  of  M-FREL  this  would  be  a  significant  damage  parameter  to  be 
modelled  as  a  pseudo-crack. 

It  was  found  possible  to  include  friction  forces  between  stringers  and 
the  sheet  but  the  accounting  was  very  complicated  with  no  practical 
improvement  in  intensity  prediction. 

The  considerable  experimental  work  in  this  project  was  divided  into 
two  parts: 

1)  Matching  da/dn  and  zK  predictions  under  constant  amplitude.  This 
led  to  flexibility  corrections  to  the  Douglas-McDonnell  formula  (121. 

2)  Further  testing  for  correlations  between  constant  amplitude  and 
random  load  data.  This  involved  best-fit  Wheeler  parameters 
associated  with  best  Forman  parameters  for  the  basic  da/dn  data. 

It  was  found  that  the  scatter  factor  of  two  recommended  by  the 
damage  tolerance  report  MIL-83444  was  not  always  conservative. 

2.2  Materials  and  Structures  Department,  Royal  Aircraft  Establishment, 
Farnborough  (Thursday  30th  April,  1987) 

I  was  met  by  Peter  Adams,  head  of  Airworthiness  Division,  MSI.  He 
said  that  Dorothy  Holford  whom  I  had  hoped  to  see  was  at  the  Madrid  meeting 
of  AGARD  (like  Salvetti)  in  particular  to  hear  a  presentation  by  a  Dr  Stuck 
from  FRG  about  characteristic  manoeuvre  sequences.  Peter  said  that  the 
method  was  oversimplified  and  that  his  feeling  was  shared  by  those  in 
Messerschmidt-Bolkau-Blohm.  I  spoke  to  Brian  Perrett  and  (Peter  Adams) 
about  the  range-pair  paper  (61  for  IFTP  and  left  some  copies.  The  number  of 
staff  at  RAE  has  fallen  and  competition  for  graduates  from  computing  firms 
has  made  recruiting  difficult. 

In  the  afternoon  I  saw  David  Rooke  who  spoke  of  intensity  solutions 
and  the  RAE  investigation  of  fretting.  To  answer  a  question  he  said  that  he 
did  not  plan  to  update  Rooke  and  Cartwright  113)  since  more  recent 
references,  such  as  Murakami  (141  were  more  comprehensive.  However  "If 
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you  read  all  the  literature  you  wouldn’t  have  time  to  write  any  books”  and  he 
thought  that  not  all  known  solutions  were  in  the  standard  compendiums.  His 
theory  on  compounding,  from  Southampton,  has  recently  been  published  as  a 
book  1151  and  1  also  saw  the  latest  summary  of  the  compounding  method  1161. 
He  also  said  the  boundary  element  method  was  attractive  and  useful  but  its 
reputation  has  suffered  by  the  propagation  of  poor  programs. 

Roger  Fenner  at  the  University  of  Southampton  is  apparently  using 
BE  methods  with  Beukner  singular  fields  117,181  and  also  Rooke,  Cartwright 
and  Aliabadi  1191. 


Dr  Rooke  promised  to  send  a  report  1201  of  the  most  recent  RAE 
results  on  fretting  fatigue  which  are  now  being  applied  to  lugs.  It  was  stated 


that  the  general  phenomenon  could  be  appreciated  in  terms  of  a  crack  in 
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In  the  notation  of  the  figure  there  are  direct  and  frictional 
components  of  intensity  proportional  respectively  to  Va  and  1/va  so  that  the 
dominant  term  depends  on  crack  length.  There  is  also  a  phase  difference 
between  frictional  and  direct  intensities  with  forces  selected  as  in 
Figure  (lc).  This  may  mean  that  additional  turning  points  are  introduced.  As 
I  mentioned,  this  can  happen  in  other  fatigue  situations  notably  the  buckling 
cracks  in  the  wing  leading  edge  of  the  CT-4.  The  report  to  be  sent  describes 
these  effects  in  detail  together  with  the  3-D  application  of  singularity 
subtraction,  related  to  a  procedure  of  L.S.D.  Morley. 

The  last  visits  of  the  day  were  to  Dr’s  Graham  Dorey  and  Chris  Peel 
of  Materials  Division  (MS4)  who  specialise  in  Composite  Materials  and 
Aluminium  Alloys.  The  composites  section  is  developing  fractography  for 
failed  composite  structures  and  working  with  Rhys  Jones  of  ARL,  using  some 
locally  developed  programs.  Dr  Dorey  mentioned  reverse  interactions  which 
could  be  tailored  to  reduce  loads  in  helicopter  blades  and  forward  swept 
wings.  Helicopter  blades  commonly  use  GFRP  rather  than  carbon  fibre  for 
superior  damage  tolerance.  Dr  Peel  is  also  interested  in  fractography  and  I 
was  also  shown  examples  of  Aluminum-Lithium  Alloys,  superplastic  forming 
with  fine-grained  material  and  precision  die  forming  (PDF).  (There  was  also 
some  superplastic  forming  at  Rockwell  International,  Section  2.7).  These 
processes  are  applicable  to  Aluminium-Lithium  Alloys  which  are  structurally 
similar  to  copper-based  alloys  with  similar  heat  treatment.  ALCAN  has  been 
licensed  to  manufacture  one  of  the  RAE  alloys  which  is  selling  well  in  North 
America  and  Europe  despite  a  price  around  21fe  times  that  of  the  alloys  being 
replaced.  Dr  Peel  said  that  the  cause  was  lack  of  competition. 

2.3  IFIP  Conference,  University  of  Aalborg  (6th  -  8th  May  1987) 

The  program  and  a  list  of  attendees  are  shown  in  Tables  I  and  n.  As 
a  gathering  by  invitation,  though  not  all  arrived,  the  IFIP  meeting  either 
included  or  mentioned  all  experts  known  to  me  and  many  more  that  were 
not.  All  the  European  oil  producers  were  represented  and  of  course  North 
America  and  Japan.  Before  discussing  the  conference  a  brief  review  1211  is  in 
order. 


In  the  modem  approach  to  structural  reliability  complete  evaluation 
of  the  risk  over  many  dimensions  is  not  normally  attempted.  It  is  assumed 
instead  that  all  variables  of  interest  (parameter  space  X  )  are  normally 
distributed,  if  necessary  after  transformation.  The  parameter  space  is  then 
dichotomised  into  safe  and  unsafe  regions  G<X}  s  0  or  otherwise  where  G  is 
usually  defined  piecewise  according  to  the  various  failure  mechanisms  to  be 
considered  (Yield  surfaces  are  typical  cases). 

It  is  then  assumed  that  one  failure  mechanism  provides  most  of  the 
risk.  Technically,  one  performs  a  Rosenblatt  transformation 

X  U  where  U  =  {U  .  .  .11} 

and  U.~  NUXO.l).  This  obviously  requires  eigenvectors  or  the 
equivalent  and  the  procedure  in  other  fields  is  far  from  new.  With  this 
mapping  the  safety  condition  becomes  g(u)  s  0;  the  reliability  index  is 
defined  as  p  =  min  (iullgtu)  «  0}  and  the  actual  risk  must  then  be  less 
that  the  normal  probability  «(-/3) . 

The  parameter  space  is  usually  a  mixture  of  geometric,  load  and 
strength  variates  and,  as  in  classical  reliability  theory,  the  mechanisms  of 
failure  are  usually  investigated  as  a  tree  of  failures,  often  combined  with 
plastic  limit  analysis. 

Current  research  and  developments  should  be  fitted  into  or  extend 
this  framework.  In  u-space  the  choice  of  B  is  effectively  a  choice  of  failure 
mode  and  open  to  engineering  judgement,  approximation  and  the  methods  of 
limit  analysis.  In  framework,  dimensions  up  to  50  are  common  so  that  it  is 
often  deemed  adequate  not  to  proceed  past  the  first  or  second  component 
failure  (FORM  or  SORM). 

Reliability,  especially  fatigue,  depends  on  time  so  that  much  effort 
is  expended  upon  the  probability  of  first  passages  of  g(U!T)  ,  say,  past  zero. 
Load  statistics,  modelling  and  estimation  were  therefore  treated  by  several 
authors  at  both  conferences. 
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One  development  that  has  excited  some  interest  is  called 
0-unzipping  1221.  In  principle  one  starts  with  any  unit  load  vector  and 
increases  it  until  first  yield  after  which  "a  giant  ten  armed  octopus" 
(Ditlevson  I2J  1  p.  1)  fixes  that  plastic  hinge  and  loading  continues.  Successive 
hinges  are  similarly  fixed  and  for  plastic  yield  lead  to  the  least  safe  ultimate 
load  vector  and  p.  At  the  conference  this  was  the  subject  of  several 
papers.  These  tried  to  extend  the  procedure,  to  brittle  failures  for  example, 
but  Ditlevson  disagreed  with  some  of  these. 

2.4  Departments  of  Civil  and  Mechanical  Engineering,  University  of 

Waterloo  (13th  -  14th  May  1987) 

At  Waterloo  most  discussions  were  with  graduate  students  and  the 
remainder  with  Professors  D.J.  Bums  and  Tim  Topper  of  the  Mechanical  and 
Civil  Departments.  I  was  told  that  their  policy  of  sandwich  instruction  has 
been  instrumental  in  attracting  students  with  the  eventual  effect  of  better 
students,  more  local  industry  and  an  increase  in  population  which  had  joined 
four  originally  small  towns.  The  main  speciality  at  the  University  was 
computer  science  (WAT  IV,  WAT  V  etc). 

David  Bums  first  showed  fatigue  tests  on  tube  joints  for  off-shore 
rigs.  Major  diameters  were  typically  1  metre  and  various  levels  of  simulation 
were  used  up  to  sea  water  environment  for  an  actual  joint.  This  was  largely 
co-operative  work  with  Denmark  (Neils  Lind  is  from  Bornholm).  Along  the 
fillet-weld  crack  growth  at  1  cm  intervals  was  measured  by  a  pulsed  DC 
potential  drop  (DCPD)  technique.  The  equipment  is  calibrated  by 
measurements  of  current  in  foil  and  boundary  element  solution  of  the  Laplace 
equation  for  current.  Conformal  mapping  was  also  used.  These  were  basically 
two-dimensional  techniques  but  the  student  describing  this  mentioned 
allowances  for  side  leakage  and  interactions. 

On  Wednesday  afternoon  I  lectured  about  the  development  of  M- 
FREL.  I  found  myself  concentrating  upon  the  Weibull  paper  I4l  with 
computing  references  confined  to  data  structure,  fitting  and  integration  of 
the  crack-damage  equations.  On  other  occasions  I  spoke  of  range-pairs  and 
the  linked-list  algorithm  for  reconstitution  15). 
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On  Thursday  I  was  introduced  to  the  Ph.D  students  Paul,  Jacques  and 
David  DuQuesnay.  The  first  project  is  the  prediction  of  life  in  45°  buttress 
threads.  From  the  time  of  Sopwith  much  has  been  known  about  local  stresses 
without  cracks  and  the  load  distribution  among  the  threads.  Paul’s  project  is 
tc  use  weight  or  Green’s  functions  (supplied  by  Jacques)  to  compute  stress 
intensities  then  crack  growth  and  life.  At  present  the  pre-crack  stresses 
assumed  for  semi-elliptic  cracks  are  far  too  conservative  because  of  the  rapid 
decay  of  concentrated  stress  away  from  the  thread-roots. 


G(x,y) 
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FIGURE  2  -  BURNS  WEIGHT  FUNCTION 


The  project  of  Jacques  was  to  extend  the  Bums  [231  weight  function 
(Figure  2)  to  edge  cracks.  This  is  exact  for  penny  shaped  cracks  and  good  for 
many  others  but  edge  cracks  present  problems.  There  is  now  a  program  which 
reduces  much  of  the  past  discrepancy  caused  by  poor  numerics. 

His  other  project  concerns  unequal  cracks  in  annuli  eg.  gun  barrels. 
Here,  if  one  crack  exceeds  the  others  these  are  shielded  but  intensity  opposite 
the  major  crack  is  increased.  When  the  cracks  are  dense  the  effect  resembles 
a  local  reduction  in  thickness  (and  effective  length)  with  a  similar  effect  at 
180°. 
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I  finally  saw  David  du  Quesnay  who  has  been  improving  a  Ford 
(Dearborn)  program  for  Neuber  analysis.  We  discussed  randomness  and  the 
Bastenaire  interpretation  of  damage  and  1  promised  to  send  him  some 
references  [24,  251. 

2.5  Institute  of  Theoretical  and  Applied  Mechanics,  University  of  Illinois, 

Urbana  (Monday  18  May,  1987) 

At  the  Talbot  Laboratory  I  was  met  by  Professor  R.T.  Shield  who  had 
arranged  visits  to  T.G.  Shawki,  F.A.  Leckie,  Kurath,  D.  Socie,  Doug  Marriot 
and  very  briefly.  Dr  K.S.  Kim. 

From  Brown  University,  Dr  Shawki  was  interested  in  plastic 
instability  in  terms  of  the  stability  of  the  governing  differential  equations  for 
plasticity.  He  said  that  the  commonly  used  linearisation  of  these  can  mislead 
about  stability  and  quoted  the  example 

Du  =  1-1  e2t|u 
!  0  -l! 

solved  by  one  Rugenius. 

Unlike  the  case  of  elastic  instability  the  reference  condition  for  the 
perturbation  equations  alters  with  time.  He  thought  that  testing  an  energy 
norm 

H>  d/dt  (u.u)  was  a  better  procedure  and  said  that  the  Lyapunov  method 
could  also  be  misleading  for  non-autonomous  equations. 

Dr’s  Socie  and  Kurath  work  together  on  multiaxial  fatigue.  For 
infinite  life  the  rules  established  by  Gough  [261  still  hold  but  finite  life  is  now 
an  area  of  research.  Experiments  at  the  Talbot  Laboratory'  use  specimens 
with  combined  tension  and  torsion  -  some  with  pressure  as  well.  Darrell  Socie 
said  later  that  their  project  was  going  well.  Prediction  was  a  matter  of 
considering  shear  and  tensile  damage  on  all  planes  in  order  to  define  that  of 
the  initial  crack  which  is  important  for  subsequent  cracking.  Miner  type 
damage  leads  to  reasonable  results  for  random  loads. 
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Professor  Fred  Leckie  began  with  an  enquiry  about  load  histories  of 
some  Victorian  rail  bridges  being  investigated  by  Paul  Grundy  at  Monash.  He 
is  a  material  scientist,  manager  of  his  department  and  has  a  student  studying 
Inconel.  His  main  interest  is  dislocation  pile-ups  and  local  stress-strain 
relations  within  grains  from  which  he  hopes  to  predict  cyclic  relations  for  the 
macroscopic  material.  He  knows  of  Brian  Cox  127,281  Section  2.7,  whose 
interests  are  similar. 

On  general  fatigue  he  said  that  "most  fatigue  can  be  designed  out", 
but  stressed  the  importance  of  statistics  and  on-site  evaluation  though  he 
acknowledged  that  this  was  impossible  for  many  mechanical  components  such 
as  shafts.  If  failures  in  these  could  ret  be  removed  by  modification  the 
Basquin  Law,  ( Ac )  1/10N  =  Const,  indicated  the  dramatic  improvements 
possible  with  reduced  loads. 

Professor  Doug  Marriot  was  interested  in  reliability  theory  or  rather 
its  failures;  this  topic  was  discussed  later  at  ICASP5  129,301.  He  Baid  that  the 
Warner  diagram  (comparing  densities  of  strength  and  load)  often  led  to  the 
false  assumption  that  these  were  always  independent  whereas,  in  fracture  for 
example,  the  metalurgical  structure  which  caused  a  low  K)c  meant  that 
corresponding  K  values  were  as  harmful  as  larger  ones  in  better  material. 
Also  in  poor  material  the  very  density  of  incipient  cracks  was  greater, 
heightening  this  effect.  He  was  also  interested  in  rogue  events  and  a  typical 
cause  about  which  he  has  written.  When  I  spoke  of  the  problems  of 
calculating  without  a  model  as  I  did  at  Tacoma  Narrows  to  Dr  Pidgeon  he  said 
that  reliability  could  still  be  estimated  although  it  became  harder.  His 
references  [31,32]  on  this  were  similar  to  the  Melchers  models  apart  from 
detailed  events.  Near  this  time  I  admitted  to  independent  damages  in  M- 
FREL  pending  the  development  of  all  the  program. 

Another  interest  of  Professor  Marriot  is  simplification  of  fracture 
mechanics  about  which  he  presents  short  courses.  Some  typical  items  for 
these  are  shown  in  Fig.  3.  Because  all  practical  cracks  have  curved  fronts 
(thumb  nail)  the  edge  connection  factor  1.1244  is  always  excessive.  Similarly 
if  a/b  for  semi-elliptic  cracks  is  finite  and  the  depth  a  <  0.4T  then  back-wall 
effects  cease  to  exist.  In  figure  (3b)  the  approximate  intensity  reflects  the 
greater  importance  of  nominal  stress  at  crack  tips. 
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I  finally  spoke  briefly  to  Dr  Kim  who  works  for  Leckie.  He  is  a 
ceramicist  who  studied  at  Brown  under  Serashi  and  he  is  engaged  in  developing 
joints  between  stabilised  ceramics  and  Inconel.  He  seemed  more  interested  in 
my  area  of  fatigue  and  we  spoke  of  polymer  fatigue,  Wolfgang  Knauss, 
fracture  mechanics  [331  and  the  reversal  in  polymers  of  temperature  and  rate 
dependence  effects  as  temperature  rises.  He  also  explained  the  role  of  cracks 
in  the  unstable  combustion  of  propellant  grains. 


2.6  School  of  Aeronautics  and  Astronautics,  Purdue  University,  West 

Lafayette,  Indiana  (Wednesday  20  May  1987) 

I  was  hosted  here  by  Professor  Alten  F.  (Skip)  Grandt  Jr.  who 
introduced  me  to  some  of  his  staff  at  Grissom  Hall  (after  Gus  Grissom)  and 
showed  some  of  the  fatigue  testing  near  Lafayette  airfield.  All  this  was  in 
the  group  Theoretical  and  Applied  Fracture  Mechanics. 

Despite  its  name  Purdue,  like  Illinois,  is  one  of  the  large  American 
state,  land-grant  universities  which  has  the  third  largest  U.S.  aeronautical 
engineering  school.  They  are  especially  proud  of  their  connection  with  the 
space  program  and  their  graduate  students  normally  include  one  or  two  RAAF 
officers.  One  of  these  is  Flight  Lieutenant  Kevin  Walker  who  proceeds  next 
to  Wamer-Robbins  AFB,  GA.  His  project  here  is  to  test  the  value  of  placard 
limitations  for  fleet  fatigue  management  using  CRACKGRO  and  experiment. 
We  also  spoke  of  the  Janus  glider  test  at  RMIT. 

I  did  not  meet  the  other  RAAF  student  Flight  Lieutenant  Adrian 
Morrison  who  is  developing  a  program  FESENT  for  computer  control  of 
testing. 


I  spoke  next  to  Dr  Tom  Farris  who  is  investigating  the  change  from 
benign  flaking  cracks  in  steel  rails  to  serious  ones  perpendicular  to  the  axis. 

The  mathematician  in  the  faculty  was  Dr  Martin  Ostoja-Starzewski 
1341  who  is  applying  random  graph  theory  to  modelling  grain  structures.  This 
resembles  the  computer  simulation  procedures  of  Cox  et  al  [27,281  whom  I 
mentioned.  Some  theory  for  this  appears  in  Vanmarcke  [351. 
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FIGURE  3  -  APPROXIMATIONS  FOR  STRESS  INTENSITY 

Until  the  arrival  of  Jack  Bogdanoff,  now  Professor  Emeritus,  I 
described  the  range-pair  paper  (6).  Martin  indicated  a  related  reference  to  be 
published  by  Macha  (361.  Bogdanoff  spoke  of  size  effects,  damage  and  related 
matters.  The  classical  Weibull  model  of  size  effect  is  a  weakest  link  theory 
of  IID  elements  in  series  which  begs  the  questions  of  defining  the  ultimate 
element  and  the  joints  between  them.  In  his  thesis  at  North  Carolina,  Pikiata 
(371  has  developed  another  model  in  which  the  weakest  flaw  is  chosen  from  a 
spatially  Poisson  density  over  the  object  of  interest.  This  has  been  applied  to 
long,  thin  continuous  fibres  and  the  length /strength  relation  for  yam  is  better 
represented  thereby. 
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This  led  to  the  topic  of  accelerated  testing  and  the  need  for  at  least 
three  levels  of  acceleration.  I  was  also  told  of  a  summary  of  fatigue  crack 
theories  to  be  presented  by  F.  Kozin  at  a  coming  SMIRT  Conference.  There 
was  finally  discussion  of  irreversible  thermodynamics  -  harder  to  follow  than 
the  size  effect  models.  The  essential  point  was  that  crack  length  was  not  the 
only  measure  of  "damage"  but  other  observables  were  also  important; 
plasticity  has  also  been  followed  and  temperature  investigated,  the  latter  by 
Luccia  and  Volta  at  ISPRA. 

At  the  airfield  laboratory  Skip  and  his  student  Ed  Tritsch  showed 
fatigue  tests  on  plexi-glass  (perspex)  models  in  which  the  visible  cracks 
allowed  photography  of  interference  fringes  for  analysis  of  crack  opening  and 
closure.  From  FEM  analysis  he  has  developed  the  formulae  for  interacting 
and  joining  thumbnail  cracks  which  drew  me  to  Purdue.  This  will  be  extended 
to  generalise  Newman's  results  1401  to  three  dimensional  or  thick-plate 
cases.  The  experimental  results  were  quite  well  predicted  from  the 
postulated  intensities  by  a  program  written  for  the  test  configuration  with 
crack  lengths  described  by  depth  and  surface  "lengths".  This  experience  here 
has  allowed  Skip  to  find  practical  estimates  of  intensity  from  nominal  stress 
distributions  and  weight  functions  and  thus  avoid  FEM  analysis  and  grid 
generation. 

Before  leaving  I  was  also  shown  a  crack  growth  test  of  grooved 
specimens  1ASTM1  and  some  results  for  Al-Li  alloy  reported  in  Fatigue  87  1411. 

2.7  Science  Center,  Rockwell  International  Corporation,  (Friday  22  May 

1987). 


Owing  to  delays  there  was  less  time  to  spend  at  Rockwell,  but 
nevertheless  I  was  met  by  Dr  Brian  Cox  in  the  morning  and  delivered  another 
lecture  about  M-FREL  which  was  well  received. 

The  Science  Center  at  Thousand  Oaks  on  US  101  is  a  separate  branch 
of  the  Rockwell  corporation  which  is  expected  to  be  profitable  in  the 
commercial  sense  though  much  of  their  income  comes  from  government 
research  contracts.  However,  much  of  their  work  is  basic,  a  counterexample 
to  those  who  would  (sic)  "make  science  relevant". 
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I  was  met  by  Dr  Brian  Cox  and  after  the  lecture  spoke  to  A.  Mike 
Mitchell,  Roberto  Kuegel,  W.  (Fred)  L.  Morris  and  John  Richardson.  Dr  Cox 
spoke  of  the  team’s  simulation  work  for  small  cracks  in  structured  materials 
and  presented  some  references  [27,281.  Mike  Mitchell  showed  me  one  of  the 
test  laboratories.  One  of  the  tests  was  in  liquid  nitrogen,  another  concerned 
metal  matrix  composites  whilst  there  was  also  automated  collection  of  da/dn 
aK  data  from  compliance  and  load  control  equipment. 


L 
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Fred  Morris  spoke  of  debonding  around  nearby  fibres  in  composites 
[431.  John  Richardson  and  I  discussed  the  life  distribution  implications  of 
crack  joining,  especially  the  relation  to  runaway  cracking  or  the  ARL-NERF 
term  rp,  continuity  of  probability  and  the  possibility  that  rapid  coalescence 
would  smooth  the  local  density  of  life.  The  context  for  this  at  Rockwell  is  the 
joining  of  many  cracks  in  heterogeneous  grains.  This  will  be  an  extension  of 
the  teams  present  work  [27,281. 


2.8  ICASP5  Conference,  University  of  British  Columbia,  Vancouver,  BC 

(25-29  May  1987) 

The  lists  of  papers  and  attendees  are  shown  in  Tables  III  and  IV. 
Because  of  the  large  numbers  the  conference  was  divided  into  three, 
sometimes  four,  separate  sessions  one  of  which  was  concerned  with  structural 
systems  as  at  Aalborg.  Other  sessions  dealt  with  geomechanics  fatigue  and 
random  dynamics.  Two  topics  with  several  papers  were  wooden  structures  and 
human  fallibility  [29,301.  In  these  ’soft’  studies  the  papers  of  Rosenblueth,  and 
Ingles  were  also  interesting.  The  first  [441  was  largely  about  the  ultimate  cost 
of  the  Mexico  City  earthquake  whilst  Owen  Ingles  1421  included  a  plea  for 
consistent  language  in  the  discipline.  Other  papers  of  special  interest  to  me 
were  Winterstein  [451  and  the  elegant  probability  theory  of  Thayaparan  and 
Karen  Chou  [461. 

In  the  area  of  fatigue,  Rackwitz  [471  and  co-workers  are 
investigating  random  processes  which  exceed  several  bounds,  related  possibly 
to  SORM. 


3.  STRUCTURAL  SYSTEMS  AND  AIRWORTHINESS 

During  the  past  decade  structural  reliability  has  developed 
considerably  and  it  now  has  its  own  jargon  (Ref.  Owen  Ingles  at  ICASP5 
[421  and  a  quite  standard  notation.  As  stated  by  Neils  Lind  at  1F1P  the 
discipline  "is  a  healthy  adolescent”  centred  on  civil  engineering  especially  as 
an  extension  of  plastic  limit  design.  The  other  technical  prop  is  the  theory  of 
continuous  random  processes  and  the  associated  statistics.  Because  of  this 
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extensive  development  and  because  the  core  of  standard  theory  assumes 
multidimensional  "structural  systems"  it  is  suggested  that  up-to-date 
airworthiness  engineers  should  adopt  the  notation  and  language  and  adapt  the 
theory. 


Some  shortcomings  remain.  These  are: 

1)  Genuinely  representative  load  cases.  Probabilities  can  be 
conditioned  for  unlikely  premises. 

2)  This  can  be  accentuated  with  &  -unzipping  and  design  code  cases  can 
be  unsuitable. 

3)  Like  mechanical  or  electronic  reliability  it  is  essentially  based  upon 
discrete  components  even  when  fatigue  is  considered.  In  structures 
this  is  typefied  by  applications  to  frameworks.  M-FREL  avoids  this 
by  considering  critical  points  or  hot-spots  rather  than  components. 

4)  Standard  theory  is  also  essentially  static.  Even  in  civil  engineering, 
brittle  failures  seem  to  be  a  problem.  Fatigued  structures  are 
treated  as  collages  from  a  predetermined  set  of  component  histories; 
again,  dynamic  changes  are  basic  to  M-FREL  14,51. 

These  matters  were  addressed  of  course,  especially  at  IFIP.  In 
fatigue  and  expecially  crack  growth  aeronautical  experience  may  be  useful  to 
structural  reliability.  The  civil  engineers  I  met  are  aware  of  this  but  some 
felt  that  airworthiness  engineers  were  mesmerised  by  retardation  models. 

It  may  be  possible  to  divide  the  field  into  different  schools.  If  so  one 
of  the  most  successful  has  the  axis  Munich,  Denmark,  VERTTEC  and  Waterloo 
which  have  contributed  to  the  development  of  PROBAN.  The  main  figures 
here  are  R.  Rackwitz,  P.  Thof t-Christensen  and  H.O.  Madsen  (501  of  VERITEC 
and  direct  effort  on  the  program  exceeds  at  least  five  man  years.  (For 
comparison  aeroelastic  programs  at  ARL  have  been  developed  for  at 
least  10  MY).  Much  of  the  theory  has  been  developed  by  Ove  Ditlevson  of 
Lyngby  who  stands  slightly  apart  from  this  school,  being  especially  critical 
of  /3  -unzipping.  Other  groups  in  the  field  are  formed  by  the  Institute  for  Risk 
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Research  at  Waterloo  (Neils  Lind)  and  the  Americans  such  as  Corotis,  J.N. 
Yang,  Wen  and  Vanmarcke. 

CONCLUSIONS 

Reliability  analysis  of  structures  is  most  developed  in  the  design  and 
assessment  of  civil  engineering  structures  particularly  oil  drilling  platforms 
and  response  to  earthquakes.  In  these  areas  the  state-of-the-art  FORM- 
SORM  technology  has  mainly  been  applied  as  a  statistical  extension  of 
standard  plastic  limit  analysis.  Most  of  the  leading  investigators  attended  at 
least  one  of  the  two  conferences  visited;  the  groups  centred  at  the  Technical 
University  of  Munich  and  the  Universities  of  Lyngby,  Aalborg  and  Waterloo 
and  Trondheim  are  important  because  they  co-operate  and  are  closely 
c  jnnected  with  VERJTEC  and  fatigue  testing  of  platform  designs.  The 
associated  names  are  shown  in  References.  The  computer  programs, 
especially  PROBAN  from  VERITEC  may  inc’ude  fatigue  but  this  seems  to  be 
done  on  the  basis  of  discrete  components  reacting  at  fixed  times  in 
accordance  with  predetermined  histories.  Therefore  there  is  a  place 
especially  in  aeronautical  applications  for  the  type  of  random  dynamic 
interactions  embodied  in  M-FREL.  However  the  practical  computation  of 
failure  risks  in  several  dimensions/components  is  far  more  advanced  for  large 
platforms  and  frames  and  I  now  envisage  procedures  similar  to  FORM-SORM 
for  computing  attrition  in  M-FREL. 

Discussions  at  the  Universities  of  Pisa  and  Waterloo  have  also  made 
it  clear  that  lengths  of  specified  cracks  do  not  fully  specify  the  degradation  of 
a  structure  so  that  a  complete  description  should  also  include  "damage" 
parameters  such  as  hole  elongation,  wear  and  corrosion  if  any  of  these 
interact  with  more  standard  parameters  in  a  way  that  varies  with  crack 
history.  The  present  structure  of  M-FREL  does  allow  this  generality  but  the 
inclusion  of  extra  cycles  from  fretting  or  buckling  effects  (Section  2.2)  is 
more  difficult. 

In  the  Laboratories  visited  there  is  little  work  being  devoted  to 
crack/crack  interaction  except  for  Rooke’s  compounding  method  and  more 
experimental  results  from  Purdue;  much  of  the  theoretical  and  survey  work 
has  already  been  published.  At  present  the  compounding  method  is  the 
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favoured  theoretical  basis  in  M-FREL.  For  crack  rate  prediction  the  journey 
yielded  fresh  information  from  Ditlevson  and  from  Cox  which  is  still  to  be 
studied. 


In  general  fatigue  four  of  the  six  laboratories  are  testing  Aluminium- 
Lithium  alloys  and  fatigue  simulations  have  started  to  allow  for  mettalurgical 
structure  (IT AM  and  Rockwell). 
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